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ABSTRACT 

R-modes in neutron stars with crusts are damped by viscous friction at the crust-core 
boundary. The magnitude of this damping, evaluated by Bildsten and Ushomirsky 
(BU) under the assumption of a perfectly rigid crust, sets the maximum spin frequency 
for a neutron star spun up by accretion in a Low-Mass X-ray binary (LMXB). In this 
paper we explore the mechanical coupling between the core r-modes and the elastic 
crust, using a toy model of a constant density neutron star with a constant shear 
modulus crust. We find that, at spin frequencies in excess of « 50 Hz, the r-modes 
strongly penetrate the crust. This reduces the relative motion (slippage) between the 
crust and the core compared to the rigid crust limit. Wc therefore revise down, by 
as much as a factor of 10^ — 10"^, the damping rate computed by BU, significantly 
reducing the maximal possible spin frequency of neutron star with a solid crust. The 
dependence of the crust-core slippage on the spin frequency is complicated, and is 
very sensitive to the physical thickness of the crust. If the crust is sufficiently thick, 
the curve of the critical spin frequency for the onset of the r-mode instability becomes 
multi- valued for some temperatures; this is related to the avoided crossings between 
the r-mode and the higher-order torsional modes in the crust. The critical frequencies 
are comparable to the observed spins of neutron stars in LMXBs and millisecond 
pulsars. 



J> 1 INTRODUCTION AND BASIC ARGUMENT neutron-neutron scattering in the nuclear matter. Their vis- 

cous damping implied that, if the neutron star core temper- 
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Gravitational radiation offers an exciting possibility to set .^^^^^ y ^ ^qS ^ ^j^j^.^^ jg ^j^g ^r the rapidly accreting 

U an upper limit on the spin frequencies of accretionally neutron "^tars in LMXBs (Brown & Bildsten 1998; Brown 

5y , spun-up neutron stars in close binaries (Bildsten 1998). 2OOO), then the r-mode instability would be triggered at a 

The source for gravitational radiation may be the mass g^ellar spin frequency just above 100 Hz. This is inconsis- 

quadrupole supported by static stresses within the neutron ^^^^ ^j^j^ observations: most neutron stars in LMXBs are 

star crust (Bildsten 1998, Ushomirsky, Cutler, & Bildsten believed to be spinning with frequencies around 300 Hz (see 

2000), or the current quadrupole from an « = m = 2 r-mode ^^^^ Klis 1999 for a review), and the existence of two 

which is driven by gravitational radiation reaction (Ander- ^ g j-adio pulsars means that rapidly rotating neutron 

sson et al. 1999, Bildsten 1998; for a review of the r-mode g^ars are somehow formed in spite of the r-mode instability 

instability see, e.g., Friedman & Lockitch 1999; for the dis- (Andersson et al. 1999). The r-mode damping rate would 

covery paper see Andersson 1998). In this paper we concen- j^^^^g be - 250 higher than that computed by Lindblom 

trate on the second option, namely that accretional spin-up (1998) to allow for such high spins 
triggers the r-mode instability, and the resulting spin-down 

torque from gravitational radiation reaction balances or ex- Lindblom et al. (1998) conjectured that, for the range of 

ceeds the accretion torque. The instability is triggered when temperatures relevant for LMXBs, mutual friction between 

the r-mode growth rate due to gravitational radiation re- ^he neutron and proton superfluids will provide the missing 

action (which is a sharply increasing function of the stellar damping. However, tour-de-force calculations of Lindblom 

spin frequency) becomes equal to the r-mode viscous damp- ^ Mendell (1999) showed that, except for ~ 2% of the ab 

ing rate lowed parameter space, mutual friction gives only a modest 

™ . 1 . ^. 1 r 1 J2 X increase of the r-mode damping rate, insufficient to account 

ine viscous damping timescale tor r-modes was first " 

i J u T- r\ A J -1 /ir>r>o\ rni fot the obscrved high spius of ueutrou stars. 

computed by Lindblom, Owen, & Morsmk (1998). They ^ ^ 

assumed that a neutron star is a fluid ball with a poly- Bildsten & Ushomirsky (2000) (hereafter BU) argued 

tropic equation of state, and that the viscosity is set by that the crust plays an important role in r-mode damping. 



© 0000 RAS 



2 Levin & Ushomirsky 



In their picture, the r-mode lives in the liquid core, and most 
of the dissipation occurs in a thin (a few cm) viscous bound- 
ary layer between the oscillating core and the motionless (in 
the rotating frame) solid crust. The damping BU computed 
is a factor of ~ 10^ higher than that of Lindblom et al. 
(1998); it implies that the r- modes would not become un- 
stable unless the spin frequency is close to break-up, which 
is significantly higher than the observed LMXB spins. An- 
dersson et al. (2000) and Rieutord (2000) have revised the 
BU estimate, down and up respectively, by factors of order 
unity. Both of these corrections are valid; the net result is 
that the BU estimate is almost unchanged. However, there is 
still a substantial correction unaccounted for in the previous 
work. In this paper, we show that, for a considerable frac- 
tion of the parameter space, the boundary layer damping 
rate is a factor of 10^ — lO'^ smaller than the BU estimate, 
and hence the critical frequencies for the onset of the r-mode 
instability are correspondingly lower, ~ 200 — 600 Hz. 

BU considered only the leading-order transverse motion 
pattern of the r-modes and assumed that these motions do 
not penetrate into the crust. They estimated the magnitude 
of the viscous coupling between these transverse motions 
and the crust, and found it to be small. Therefore, their 
picture is that of a liquid in a bucket - the amplitude of 
the relative motion (slippage) between the crust and the 
core is just equal to the amplitude of the fluid motion at 
the boundary. However, small radial motions of the r-mode, 
neglected by BU, are very effective at coupling the crust 
and the core. Above a small critical frequency, r-modes do 
strongly penetrate into the crust, and the crust and the core 
oscillate almost in unison with very similar amplitudes. The 
slippage Su in this regime is only a small fraction of the total 
r-mode motion, u, with typical values of Su/u ~ 0.05 — 0.1, 
which results in a correction of a factor of {5u/uf to the 
BU viscous boundary layer damping rate. 

The fact that the crust is not rigid is central to our 
story. In fact, the crust is more like a gel: the shear modulus 
is much smaller than the bulk modulus. The I = m = 2 
torsional mode in a non-rotating crust, which has the same 
angular displacement pattern as the corresponding r-mode, 
has a frequency /„ — 50 Hz (McDermott et al. 1988). This 
frequency is a few times lower than the r-mode frequency 
in rapidly rotating stars, indicating that the elastic restor- 
ing force is quite weak compared to the Coriolis force which 
is responsible for the r-mode oscillations. Therefore, at suf- 
ficiently high spin frequencies, one would expect the crust 
to oscillate more or less like a liquid, with elasticity only 
slightly modifying the mode structure. Indeed, if the shear 
modulus ^ of the crust were exactly zero, i.e., if the crust 
were completely fluid, there would be no slippage 5u between 
the crust and the core. Since /i is non-zero but small, one 
would expect the slippage to be proportional to the ratio of 
the elastic restoring force to the Coriolis force, i.e., roughly 
5u/u ~ (/cr//rmodc)^- However, our numerical work shows 
that the slippage is not a monotonically decreasing function 
of the spin frequency; it rises sharply at the frequencies of 
avoided crossings between the r-mode and the higher-order 
torsional modes of the crust. Figure |l| shows the slippage as 
a function of the spin frequency computed for two neutron 
star models with different crust thicknesses (10% and 20% 
of the stellar radius) . 

We begin in § |2| by outlining the equations governing the 



pulsations of constant density stars with fluid cores and solid 
crusts and presenting the numerical results for the depen- 
dence of the slippage 5u/u on the spin frequency. In § ra we 
use the computed slippages to construct the critical r-mode 
stability curves; see Figure ^. The resulting curves are quite 
sensitive to the thickness of the neutron star crust. In par- 
ticular, for a sufficiently thick crust, the stability curve has 
a peculiar multi-valued structure, which is a consequence 
of the non-monotonic dependence of the relative slippage 
on the spin frequency. The stability curve determines the 
maximum spin of a neutron star as a function of temper- 
ature; we find that it generally goes through the region of 
observed spin frequencies of neutron stars in LMXBs and in 
millisecond radio pulsars. Broad implications of our results 
are discussed in § ^. 



2 CORE-CRUST BOUNDARY SLIPPAGE 
2.1 Basic formalism 

Oscillation modes of a rotating neutron star with an elastic 
crust were analyzed by Lee & Strohmayer (1996) (hereafter 
LS). Their analysis was performed for the case of slow rota- 
tion, which was treated as a perturbation to other restoring 
forces, and the oblateness of the equilibrium configuration 
was neglected. We broadly utilize their general formalism, 
however we derive equations specialized to the constant- 
density case. For r-modes the correction to the eigenfre- 
quency due to the non-sphericity of the equilibrium con- 
figuration is of the same order as the correction due to the 
radial motion (Provost et al. 1981). However, we expect that 
the oblateness of the star will not change the qualitative pic- 
ture for the crust-core slippage: at high enough spin frequen- 
cies, the core and the crust will move almost in unison. The 
quantitative effects of non-sphericity are not treated in this 
paper. 

For our toy model, we assume that the outer part, r^i < 
r < R, of the neutron star of radius R possesses a small 
shear modulus /i(r), such that /^(r) = const. The neutron 
star fiuid is assumed to be incompressible, has a constant 
density p, and is rotating with angular frequency fl. 

The general displacement u of a fluid element, in a frame 
rotating with the neutron star, is given by (Zahn 1966, LS): 



i?=^(^i,(r)yi^f-f i;,(r)VYi^ + w,(r)f x VY,^)e"^\ (1) 
I 

where uj is the mode frequency in the rotating frame, m is 
the fixed azimuthal number, and VYim is the gradient of 
the spherical harmonic evaluated on a unit sphere. We are 
interested in the I — m r-mode. For a purely fluid star, to 
0{n/ni,f, where fik = {GM/R^f/^ is the keplerian an- 
gular velocity, the displacement for this mode is given by 



{wif X VYii + Ui+iWi+i I + Ui+iY"i+i if) e" 



(2) 



i.e., vi+i/wi and ui+i/wi are of order (Sl/f^k)^. In the solid 
crust, the displacement vector need not satisfy the ordering 
assumed in Eq. However, at least for modes with ui close 
to 2n/3, the Coriolis force is the dominant restoring force, 
and we expect vi+i/wi and ui+i/wi to be small; this ordering 
is confirmed by our numerical results (see § H). 
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We now define dimensionless frequencies Q — Q/fl]^ and 
uj = a;/f2k, as well as the dependent variables similar to those 
utilized by LS: 

(3a) 
(3b) 
(3c) 
(3d) 
(3e) 

(3f) 

In the constant density case, Z2 is not a dynamical variable, 
as it can be eliminated using the continuity equation, V-u = 
0, or 
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(4) 



However, even in the p — const case, the Eulerian pres- 
sure perturbation is non-zero, and, to leading order, is given 
by Sp — 5pi+i{r)Yi+i i. We define the dimensionless pres- 
sure perturbation, p = Spi+i/pgr, where g = GM{r)/r^ 
is the local gravitational acceleration. Finally, we define 
the dimensionless shear modulus, jl = p/pgr. Through- 
out this paper, we use fl = 1.5 X 10'"^{R/rf, chosen so 
that the the frequency of the I = m = 2 fundamental 
torsional crustal mode is close to the ~ 50 Hz value com- 
puted by (McDermott et al. 1988) for a non-rotating neu- 
tron star; the shear wave speed corresponding to this choice 
is ct = {p/pf^ = 1.7 X 10* cm 8"^ 

These variables allow us to write the equations of mo- 
tion in the crust as a system of first-order differential equa- 
tions: 



dzi 
dr 
dp 
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= {l + l){l + 2)z3-3zi 
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where jj^i = (2^ + 3)"^'''^. The derivation of these equations 
is analogous to that of Eqs. (9)-(21) of LS. In the liquid core, 
the only nontrivial equations are 

dp 



(6a) 



dz\ ,, ^, 

r-r^ = -/-I- 4)21 
dr 
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and 
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(7) 



(see, e.g., Kokkotas & Stergioulas 1999 or LS). 

Both at the crust-core boundary, and at the top of the 
crust, we require the tangential traction to be zero, i.e.. 



24 = 26 = 0. 



(8) 



We hence neglect the viscous coupling between the trans- 
verse r-mode motions and the crust. This is the coupling 
considered and found negligible by BU. At the crust-core 
boundary the radial displacement 21 is continuous, and so is 
the radial component of the traction, Ap — 2pdui+i/dr, or 

zi{r~) = 2i(r-+) (9a) 
p{rt) - m[2(Z + 1)(Z + 2)23 (r,+ )- 421] =j5(rr). (9b) 

At the top of the crust, the radial component of the traction 
is zero: 

p{R) - zi{R) - ji[2{l + 1)(Z + 2)23(7?) - 421 (i?)] = (10) 

Equations (^ — (jl^) provide the boundary conditions for the 
equations of motion. We normalize our solutions to have unit 
transverse displacement just inside the liquid core. 



ZbiVc ) = 1- 



(11) 



With this normalization, the slippage is 5u/u = [{z^ij-'^) — 
25(r^))^ + {zz{r~) — 23(7-^))^]"'"''^ (where we used our as- 
sumption, 23 <C 25). 

2.2 Numerical results 

The keplerian spin frequency fik, the shear modulus p, and 
the ratio of the crust thickness to the stellar radius rd R are 
the only parameters necessary to describe a constant density 
star (we neglect the thin surface ocean and the modes that 
are present in it, as they are unlikely to influence the crust- 
core slippage). The crust thickness is roughly proportional 
to the scale height at the crust-core boundary, and is hence 
(X R? jM. In order to survey the range of parameters, we 
carry out our calculations for two models, a "thin" crust 
model with r^jR ~ 0.9, and a "thick" crust model with 
rc/ R — 0.8. Because of the differences in the normal mode 
spectrum between these two cases, we shall see that the r- 
mode instability curves for the two models have qualitatively 
different behavior in the temperature range of interest for 
accreting neutron stars, T = 10* — 10^ K. 

The left panel of Figure |l| displays the results of our 
calculations for I — m — 2 for the thin crust model, while the 
right panel shows the results for the thick crust model. Let us 
first concentrate on the top plots, which show the dispersion 
relation, aj(f2), for the two models. Contrary to the case of a 
purely fluid star, where, in the isentropic case, there is only 
one "classical" I — m r-mode, all the modes displayed in 
Figure ^ have displacements of the form given by Eq. (^, 
with Wl (X m the core. However, only the mode with uj 
closest to the "classical" r-mode value, 20/3, is of interest to 
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Figure 1. Top left; Mode angular frequency uj (in the rotating frame) as a function of the angular spin frequency O for the five lowest- 
frequency modes of the thin crust model. The dashed line shows the "classical" r-mode dispersion relation, ui = 20/3, and the dots mark 
the mode closest to this relation, which has the smallest slippage. Bottom left: absolute value of the slippage of the mode marked with 
thick dots in the top panel (i.e., the minimum slippage) as a function of Q. Right: same as top left, but for the thick crust model. 



us. This mode is denoted by thick dots in Figure When 
Lo is much different from 2f2/3, the displacements in the 
crust are much bigger than the displacements in the core, 
resulting in large slippage, and, hence, large damping. In 
essence, unless the mode frequency is close to the r-mode 
value, the elasticity of the crust plays the dominant role in 
determining the mode properties, and the mode is a crustal 
mode, rather than a core r-mode. 

Consider the low spin frequency (f2 <^ {fi/ pR^Y''^ ~ 
COSfik) hmit of our results. In this case, the mode with 
the lowest frequency in Figure ^ is the "classical" I = m 
r-mode, which has the transverse displacement W2 oc in 
the core and negligible displacements in the crust. This be- 
havior (liquid in a bucket) is the approximation used by BU. 
The modes at higher frequencies are torsional crustal oscil- 
lations, which have appreciable displacements in the crust, 
and negligible displacements in the core. The fundamental 
(nodeless) torsional mode has lj = 0.03f2fc for both models. 
The lack of sensitivity of the frequency of this mode to the 
thickness of the crust was first noted by Hansen & Ciofh 
(1980). The modes at higher frequencies are overtones, and 
we note that the mode spectrum for the thick crust model is 
denser than that of the thin crust model. Qualitatively, the 
mode frequency is determined by fitting a certain number 
of nodes of the eigenfunction between the boundaries of the 
crust, and hence a thicker crust supports more modes. 

As the spin frequency is increased beyond fl « 
{Ijl/ pB?)^^^ , the behavior of the modes changes. The fre- 
quency of the core r-modes rises to meet with the frequencies 



of the modes which originally resided in the crust, resulting 
in avoided crossings. At the first avoided crossing, the core 
r-mode is coupled with the lowest-order torsional mode of 
the crust; at this avoided crossing the core r-mode strongly 
penetrates the crust. As the spin frequency increases beyond 
the first crossing, the crust-core slippage decreases along the 
branch closest to the "classical" r-mode frequency, u) = 2f2/3 
(marked by a dashed line in Figure |l|). Since the boundary 
layer damping rate is proportional to {5u/uf , this branch is 
most interesting for us as it corresponds to the r-mode with 
the lowest damping, and, hence, the most unstable. 

The relative slippage as a function of the spin frequency 
is plotted in the bottom panels of Figure ^. We always plot 
slippage of the r-mode from the branch with the lowest 
damping, marked by thick dots in Figure |l| (i.e., the most 
unstable r-mode). As the spin frequency increases beyond 
the first crossing, the slippage decreases to « 0.1. How- 
ever, the thick crust model has another avoided crossing 
at n ~ 0.25r2k, and the slippage rises at this frequency be- 
cause the core r-mode couples resonantly with a higher-order 
torsional crustal mode. In the thick crust model, the "val- 
ley" in which the second mode dominates only extends to 
Q, ~ 0.25r2k- After the avoided crossing, the third mode 
has the smallest slip, 4 x 10~^ < 5u/u < 1. Note that the 
slip value is ~ 1 near the avoided crossing and is ^ 0.1 
for 0.25r2k < < 0.3f2k, or for spin frequencies between 
~ 500 and 600 Hz. In this regime, the viscous boundary 
layer damping rate approaches the high BU value, and so 
the critical spin frequency for the r-mode instability is high 
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as well. As we now show, the shape of the r-mode stability 
curve depends quite sensitively on the thickness of the crust. 



3 CRITICAL STABILITY CURVES 

The dissipation rate in the viscous boundary layer is set 
by the amplitude of the relative motion between the core 
and the crust. BU assumed that the crust remains static 
in the rotating frame, and hence the amount of boundary 
slippage Su/u = 1. However, as shown in the previous sec- 
tion, the slippage Su/u is generally significantly less than 1. 
Therefore, the BU estimate of the damping rate has to be 
multiplied by {Su/u)"^ . The corrected r-mode damping rate 
is [cf. Eq. (4) of BU] 



(12) 



Tvbi MiaTs p VkHz/ 

where / = Q,/2-n is the neutron star spin frequency, R^, 
MiA and Tg are the radius, the mass and the core temper- 
ature of the neutron star in units of 10 km, IAMq, and 
10* K, respectively, p is the density at the crust-core in- 
terface, pb = 1.5 X 10^'* g cm~^ is the estimated value of 
this density (Pethick et al. 1995), and F is the fitting pa- 
rameter used by Cutler & Lindblom (1987) to parameterize 
shear viscosity found by Flowers & Itoh (1979). The Ekman 
layer correction introduced by Rieutord (2000) for the con- 
stant density case is of order unity and is not included in 
Eq. (0); see Lindblom et al. (2000) for a treatment appli- 
cable for more realistic, non-uniform density stars. We also 
neglected the contribution of the elastic energy of the crust 
to the total mode energy, as well as the small deviation of 
the mode frequency from u = 20/3. 

The r-mode growth rate due to gravitational radiation 
is given by 



— ~ 0.02 s'^MiARi (t4^) 



VkHz/ 



(13) 



Note that this expression is different from the one used by 
BU. Andersson et al. (2000) pointed out that BU incorrectly 
rescaled rgw used by Lindblom et al. (1998) to the canon- 
ical R = 10 km. Andersson et al. (2000) gave their own 
estimate of the r-mode growth rate for a constant density 
star. Eq. (^ is the correct Tgw for an n = 1 polytrope, ne- 
glecting the deviation of the transverse displacement from 
wi oc r' in the crust, as well as the small deviation of the 
mode frequency from (jj — 20/3. 

While we computed the slippage using a constant den- 
sity model, we utilize Tvbi (Eq. Q) and rg„ (Eq. [Q) ap- 
propriate for a more realistic n = 1 polytrope. Compared 
to more realistic models, using Tgw and r^rmvu for a con- 
stant density model overestimates the gravitational radia- 
tion growth rate (by putting more mass at larger radii), and 
underestimates the strength of the boundary layer damping 
(by virtue of a larger mode energy). Both of these effects re- 
sult in lower critical frequencies for the onset of the r-mode 
instability. Slippage calculations for more realistic models 
are in progress. 

The stability curve is obtained by equating the 
gravitational-radiation-induced growth rate with the viscous 
dissipation rate. Since we are primarily interested in temper- 
atures below the melting temperature of the crust, ~ 10^" K, 
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Figure 2. Critical frequencies as functions of the core temper- 
ature for the modes in the thin crust model (solid line) and the 
thick crust models (dotted line). For comparison, we also dis- 
play the critical frequencies for the n = 1 polytrope fluid model 
(dashed line), and for a model with a crust and Su/u = 1 inde- 
pendent of the spin (dot-dashed line) . 

we neglect the dissipation due to bulk viscosity. We also ig- 
nore the dissipation due to the shear in the interior of the 
star, as well as any potentially important dissipation mech- 
anisms that may be operating in the crust, and presume 
that the viscous boundary layer is the dominant energy sink. 
Since we do not have an analytical expression for the slip- 
page, it is more convenient to evaluate the critical tempera- 
ture, Tg, such that Tg 



10 



Tvbi: 



Ta = 0.5 



F'^' Pb I f 



RlMl 



p VkHz 



(14) 



Figure g shows the critical stability curves for models with 
thin (solid line) and thick (dotted line) crusts, with F, Mia, 
i?6 and pb/p all equal to one. In other words, the proper 
label for the abscissa of this plot is TRImIaF~^'^ p/ Pb- For 
comparison, we also plot the stability curve for a purely fluid 
n — 1 polytrope (Lindblom et al. 1998, Lindblom et al. 1999, 
dashed line), as well as the curve with Su/u = 1 independent 
of spin frequency (dash-dotted line, cf. BU). Neglecting the 
shear viscosity in the interior of the star is clearly justified. 

For the thin crust model (solid line), the slippage in the 
frequency range of interest does not rise or fall sharply, and 
so the critical frequency is a single- valued, monotonically de- 
creasing function of T. However, for the neutron star with 
a thick crust, the slippage and the damping rise sharply 
at the spin frequency close to the second avoided crossing. 
Therefore, the critical temperature is not a monotonically 
decreasing function of the spin frequency /, and the stabil- 
ity curve is multi-valued. Such a non-trivial stability curve 
should be interpreted in terms of the critical temperature. 
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rather than the critical frequency: for temperatures larger 
than the critical one given by Eq. (^^, i.e., for points in 
the T — / plane which lie to the right of the viscous bound- 
ary layer part of stability curve, the r-mode is unstable and 
grows. Otherwise, the r-mode is stable and decays. 



4 DISCUSSION, BUT NO BETS 

Observationally, the spin frequencies of neutron stars in 
LMXBs have been inferred from the QPO-like flux mod- 
ulation during the type I X-ray bursts (the so-called burst 
QPOs, see van der Klis 1999 for a review). In two out of 
the six sources displaying the burst QPOs, the frequency 
is close to 300 Hz, while in the remaining four it is close to 
600 Hz. In one of the latter four systems. Miller (1999) found 
a subharmonic at around 300 Hz, and has argued that this 
subharmonic corresponds to the true spin frequency. Same 
relation is sometimes believed to hold for the other ~ 600 Hz 
bursters. In addition, « 300 Hz spin is often argued on the 
basis of a beat frequency interpretation of the kHz QPO 
peak separation in all six of the burst QPO systems, as well 
as for the neutron stars in thirteen additional LMXBs that 
do not display burst QPOs. On the other hand, there is no 
compelling theoretical picture for why the subharmonics of 
the ~ 600 Hz burst QPOs should be so weak or invisible, so 



a bimoa ial apin diatribution ia certainly not ruled out. 

Th e now critical stability curvoB computed in thic pa 



per for two toy models of neutron stars cut right through 
the observed range of spins both for neutron stars in LMXBs 
and for millisecond pulsars, for temperatures of interest for 
accreting neutron stars (between 10* K and 10® K; Brown 
& Bildsten 1998, Brown 2000). An accreting neutron star 
can continue its spin-up so long as it is located to the left of 
the stability curve in Figure ^ However, once it crosses the 
critical stability curve and is located to the right from it, 
the r-mode will grow and halt the spin-up, probably forcing 
the neutron star into a thermal runaway cycle (Levin 1999; 
Spruit 1999). The details of what happens during this run- 
away and how violent it is are a topic of current research. 
For example, the displacements induced in the crust may 
cause it break when the strain (which is of order the r-mode 
amplitude a, as defined in Lindblom et al. 1998) exceeds a 
critical value (< 10~^, see § 6.1 of Ushomirsky et al. 2000 
and references therein for a summary of estimates of crustal 
breaking strains) or melt if the heating due to the viscous 
boundary layer raises its temperature above ~ 10^° K (see 
Lindblom et al. 2000 for the calculation of the r-mode ampli- 
tude required to melt the crust). The evolution of the spin 
frequency and the temperature depends on the saturation 
of the r-mode growth, which is not currently understood. 
Despite these uncertainties, it is clear that, if the r-mode in- 
stability operates in accreting neutron stars, the non-trivial 
shape of the stability curve and its sensitivity to the crustal 
structure will be reflected in a peculiar, non-uniform spin 
distribution of accreting neutron stars and millisecond pul- 
sars. 

Much work is necessary in order to make robust compar- 
isons with observations. R-mode calculations for more realis- 
tic neutron star models, taking into account the complicated 
structure of the crust are in progress. We have not consid- 
ered the possible superfluid nature of neutrons and protons 



at the crust-core interface. In particular, we have not consid- 
ered the pinning of the superfluid vortices to the crust. We 
assumed that the boundary layer is laminar, which may not 
be true for large r-mode amplitudes. Signiflcant magnetic 
field at the crust-core boundary may also affect the crust- 
core coupling (RezzoUa et al. 2000; Ho & Lai 2000). Thus, 
our estimate of the r-mode damping may change when more 
complex physics is considered. 
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